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^ ' Abstract: We study a matrix model with a cubic term, which incorporates both the 

' fuzzy X and the fuzzy as classical solutions. Both of the solutions decay into the 

vacuum of the pure Yang-Mills model (even in the large-A^ limit) when the coefficient of 
the cubic term is smaller than a critical value, but the large- behavior of the critical point 
is different for the two solutions. The results above the critical point are nicely reproduced 
by the all order calculations in perturbation theory. By comparing the free energy, we find 
that the true vacuum is given either by the fuzzy or by the "pure Yang-Mills vacuum" 
depending on the coupling constant. In Monte Carlo simulation we do observe a decay of 
the fuzzy x into the fuzzy at moderate N , but the decay probability seems to 
be suppressed at large N . The above results, together with our previous results for the 
fuzzy CP^, reveal certain universality in the large- A'^ dynamics of four-dimensional fuzzy 
manifolds realized in a matrix model with a cubic term. 
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1. Introduction and summary 

Fuzzy spheres [1], or fuzzy manifolds in general, have recently attracted much attention 
in various branches of particle physics. One of the motivations comes from the general 
expectation that the noncommutative geometry, which is characteristic to those manifolds, 
provides crucial links to string theory and quantum gravity. Indeed Yang-Mills theories 
on noncommutative geometry appear in a certain low energy limit of string theory [2]. 
There is also an independent observation that the space-time uncertainty relation, which is 
naturally realized by noncommutative geometry, can be derived from some general assump- 
tions on the underlying theory of quantum gravity [3]. Another motivation is to use fuzzy 
manifolds as a regularization scheme alternative to the lattice rcgularization [4]. Unlike 
the lattice, fuzzy manifolds typically preserve the continuous symmetries of the space-time 
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considered, and hence it is expected that the situation concerning chiral symmetry [5-19] 
and supersymmetry might be improved. 

As expected from the connection to string theory [20] , fuzzy manifolds appear as clas- 
sical solutions in matrix models with a Chern-Simons-like term [21-25] and their dynamical 
properties have been studied in refs. [26-39] . One can actually use matrix models to define 
a regularized field theory on the fuzzy spheres as well as on a noncommutative torus [40], 
which enables nonpcrturbative studies of such theories from first principles [41]. These ma- 
trix models belong to the class of the so-called dimensionally reduced models (or large- 
reduced models), which is widely believed to provide a constructive definition of superstring 
and M theories [42,43]. In fact there are certain evidences in the JIB matrix model [43] 
that four- dimensional space-time is generated dynamically [44 47]. In rcfs. [44-46] the 
free energy of space-time with various dimensionality has been calculated using the gaus- 
sian expansion method, and the free energy turned out to take the minimum value for 
the four-dimensional space-time. In ref. [47] it was found that the fuzzy x (but not 
the fuzzy S^) is a solution to the 2-loop effective action. See also refs. [48-60] for related 
works on this issue. The fuzzy sphere is also useful [61] in the Coset Space Dimensional 
Reduction [62,63]. 

In view of the variety of contexts in which fuzzy manifolds have been studied, we 
consider it important to study their nonpcrturbative dynamics from first principles by 
Monte Carlo simulations. In ref. [64] we have studied the dimensionally reduced 3d Yang- 
Mills-Chern-Simons model, which incorporates the fuzzy S^ as a classical solution [22] . We 
have observed a first-order phase transition as we vary the coefficient (a) of the Chern- 
Simons term. For small a the large- behavior of the model is the same as in the pure Yang- 
Mills model, whereas for large a a single fuzzy appears dynamically. The emergence of a 
fuzzy sphere in matrix models may be regarded as a prototype of the dynamical generation 
of space-time since it has lower dimensionality than the original dimensionality that the 
model can actually describe. In ref. [65] we have performed the all order calculations in 
perturbation theory around the fuzzy S^ solution following the proposal in ref. [34], and 
confirmed that the Monte Carlo results for various observablcs in the "fuzzy sphere phase" 
can be nicely reproduced. For obvious reasons it is interesting to extend these works to 
/our-dimensional fuzzy manifolds. In ref. [66] we have studied a matrix model incorporating 
the fuzzy S^, and in ref. [67] a matrix model, which incorporates the fuzzy CP^ as well as 
the fuzzy S^. 

In this paper we study the fuzzy S^ x S^, which has been studied extensively in the 
literature [32,34,36-39,47]. We perform perturbative and nonpcrturbative studies of a 
6d Yang-Mills model with a cubic term, which incorporates the fuzzy S^ x S^ as well as 
the fuzzy S^ as classical solutions. Both of the solutions become unstable (even in the 
large- AT limit) when the coefficient of the cubic term is smaller than a critical value, but 
the large-AT behavior of the critical point is different for the two solutions. The results 
above the critical point are nicely reproduced by the all order calculations in perturbation 
theory. By comparing the free energy, we find that the true vacuum is given either by the 
fuzzy or by the pure Yang-Mills vacuum depending on the coupling constant. In Monte 
Carlo simulation we do observe a decay of the fuzzy S^ x S^ into the fuzzy S^ at moderate 
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N, but the decay probability seems to be suppressed at large A''. 

In fact the above results are qualitatively similar to the fuzzy CP^ case as far as large- 
N properties are concerned. This reveals certain universality in four-dimensional fuzzy 
manifolds realized in a matrix model with a cubic term. In the case of the fuzzy [66] 
we had to add a quintic term, which led to a totally different situation. We note, however, 
that even within matrix models with a cubic term, we have found two directions which lead 
to qualitatively different results. In ref. [68] it is shown that an additional "mass term" 
can make various fuzzy sphere solutions almost degenerate at the classical level, and it is 
possible that the true quantum vacuum is described by a set of coincident fuzzy spheres, 
giving rise to nontrivial gauge groups. In ref. [69] we have shown that super symmetry 
removes the quantum effects, and as a result the single fuzzy sphere becomes always stable 
if the large- A?^ limit is taken in such a way that various correlation functions scale. 

This paper is organized as follows. In section 2 we define the model and discuss its 
classical solutions. In sections 3 and 4 we study the properties of the fuzzy x and the 
fuzzy S^, respectively, by performing Monte Carlo simulations and the all order calculations 
in perturbation theory. In section 5 we determine the true vacuum of the model based on 
the comparison of the free energy. In section 6 we discuss a transition from the fuzzy 
X to the fuzzy observed in Monte Carlo simulation. In the appendices we provide 
the details of our calculations. 

Note added: Part of this work has been reported by one of the authors at YITP 
workshop "Quantum Field Theory 2004" , July, 2004, at a meeting of the Physical Society 
of Japan, Kochi, Sep. 2004 and at Dublin Institute for Advanced Studies, Dec. 2004. 
While we were preparing this article, we received a preprint [39], which discusses the phase 
diagram of a similar model based on the one-loop effective action. 

2. The model and its classical solutions 

The model we study in this paper is defined by the action 



where [ji = 1, • • • ,6) are traceless N x N hermitian matrices. Here and henceforth 
we assume that repeated Greek indices are summed over all possible integers. The rank-3 
tensor f^^p in the cubic term is given by 
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for /X, z^, p = 1, 2, 3 , 
for fi,^, p = 4:,5,6 , 
otherwise , 



(2.2) 



where Cf^^p is a totally antisymmetric tensor with 6123 = €453 = 1. 
The classical equation of motion reads 
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']] -i^Uvp [Ay^Ap] = . 



(2.3) 
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In addition to the commutative solution, which exists also for a = 0, it has the fuzzy x 
solution, which is defined by 



A., 



-v^(ni,n2;fc) 



dcf 



a Ir 



L 



Ife for /x= 1,2,3, 
Ife for /X = 4, 5, 6 , 



(2.4) 



where L^"^ represents the n-dimensional representation of the SU(2) Lie algebra, and the 
integers ni, 77-2 and k should satisfy N = nin2 k. Let us define the "Casimir operators" 
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1 = : Q2 



(2.5) 



Plugging the solution (2.4), we obtain Qj = RjIn, where Rj = \a-\/ (nj)^ — 1 for j = 1, 2. 
This classical solution (2.4) therefore describes k coincident fuzzy x with the radii 
R\ and R2 in the 123- and 456- directions, respectively. If wc expand the model around 
the solution, we obtain the U(A;) gauge theory on the fuzzy X generalizing the work 
on the fuzzy [22]. The classical action for the solution is given by 



(2.6) 



In what follows we focus on the symmetric fuzzy x (n\ = 77-2) and the fuzzy (77-2 = 1). 
Monte Carlo simulation is performed using the heat bath algorithm as in refs. [49,64]. 
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Figure 1: The obscrvablcs obtained by Monte Carlo simulation with the fuzzy x start are 
plotted against a = aN* for N — 16,25,36 (i.e., n = 4,5,6). The dotted (dashed) lines represent 
the classical (one-loop) results at large A''. The solid lines represent the all order results at large N 
obtained above the critical point. 



3. Properties of the fuzzy x 

In this section we study the properties of the single {k = 1) fuzzy x S^. We perform 
Monte Carlo simulation taking = x^"'"'^'* as the initial configuration, where n = y/N . 
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Let us consider the expectation value of the action (2.1) and the "extent of space-time" 
tr(^^)^ = tr {Qi + Q2)- In figure 1 we plot these quantities obtained by the simulation 
with = 16,25,36 against a = aN^. We observe a discontinuity at a ~ 2.5, which 
implies the existence of a first-order phase transition. The critical point agrees with the 
result (B.25) obtained analytically from the effective action. 

Above the critical point we can calculate the observables in the large-iV limit to all 
orders in perturbation theory as 



1 a'^ 5 8 448 3520 

^(5)--y^ + 2 + ^ + ^ + -T^ + --- > (3-1) 

^ tA\A _ 4 _ 40 _ 768 _ 18304 _ 



(See appendix B.4 for derivation.) This result, as well as the classical and one-loop results, 
is plotted in figure 1. We observe that the Monte Carlo data for a > tier approach the all 
order results as N increases. 
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Figure 2: The observables obtained by Monte Carlo simulation with the fuzzy start are plotted 
against a = aN^ for N = 16,25,36. The dotted (dashed) lines represent the classical (one-loop) 
results at large TV. The solid lines represent the all order results at large A'' obtained above the 
critical point. 



4. Properties of the fuzzy 



Next we study the properties of the single {k = 1) fuzzy S^. We perform Monte Carlo 
simulation taking Aj^ = xj^^'^'^^ as the initial configuration. In figure 2 we plot the results 

for N = 16, 25, 36 against a = aN^ . We observe a gap at the critical point a ~ 3.0, which 
agrees with the result (C.8) obtained from the effective action. The Monte Carlo results 
above this point agree well with the all order results in perturbation theory at large N , 
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which are given by 
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(4.1) 
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(See appendix C.3 for derivation.) 



5. The true vacuum 



In this section we determine the "true vacuum" by comparing the free energy, which can 
be obtained to all orders in perturbation theory as 
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^feS2xS2 
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(5.1) 
(5.2) 



for the k coincident fuzzy x and the k coincident fuzzy S^, respectively. (See appen- 
dices B.4 and C.3 for derivation.) Note that these solutions are stable above the critical 
points given by (B.25) and (C.8). 

Let us first compare the free energy among different values of A;. In figure 3 we plot 
the free energy for the fuzzy S^ x S^ (left) and S^ (right) against a and a, respectively, for 
various k. The free energy for A; > 1 is always larger than that for A; = 1, which implies 
that the dynamically generated gauge group is U(l) for both the fuzzy x S^ and the 
fuzzy cases. 
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Figure 3: The free energy for the k coincident fuzzy S'^ x (left) and the k coincident fuzzy 
(right) is plotted for fc = 1, 2, • • • , 6. We have subtracted the irrelevant constants (4 log N and 
5 log A'', respectively) to make the quantity finite in the large- A'' limit. 

We therefore take A; = 1 and compare the free energy for the fuzzy x and the fuzzy 
S . In order for the single fuzzy S x S to be stable, we need to have a > acr ■ In 
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that regime, however, a ~ 0{N^^'^^ and therefore W^^^'sp. ~ — 0(A^^), meaning that the 
single fuzzy has much smaller free energy than the single fuzzy x S^. 

As we have done in ref. [68], we can also calculate the free energy WVm in the Yang- 
Mills phase for a in the large- AT limit since it should agree with the free energy for 
the pure Yang- Mills model [a = 0) studied in ref. [71]. We obtain (See Appendix D.) 

-^HVM = 31ogA + 6|-0.33 + log(337r-^)} . (5.3) 

Comparing this result with (5.2) for A; = 1, we find that the true vacuum of this model is 
given by the single fuzzy for a > acr and by the pure Yang- Mills vacuum for a < 
where the critical point is 

acr^(481ogAj . (5.4) 

This result is exact at large N since the calculation of VFjfc=i 32 and Wym are both reliable 
at the critical point Ocr- 

Note that the lower critical point Ocr = air~^^ \ at which the single fuzzy becomes 
unstable, is below ttcr- On the other hand, from Monte Carlo simulations we find that the 
upper critical point, at which the pure Yang-Mills vacuum becomes unstable, is ~ 1.51, 
which is above acr = ^^'^cr- The existence of such three kinds of critical points is typical 

to first-order phase transitions. In the region Ocr < a < we observe a hysteresis 
behavior in simulations similarly to the one observed in ref. [64]. 



6. A transition from the fuzzy 
X to the fuzzy 

For a > CKcr, since the true vacuum is 
given by the fuzzy S^, we may be able 
to see the fuzzy x prepared as 
the initial configuration decay into the 
fuzzy in Monte Carlo simulations if 
we make a very long run. 

In order to distinguish the classical 
solutions, it is convenient to look at the 
eigenvalues of the "Casimir operators" 
Qi and Q2- In figure 4 we plot the eigen- 
values obtained for N = 9 and a = 2.6 
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Figure 4: The eigenvalues of the "Casimir opera- 
tors" Qi and Q2 arc plotted against the number of 
the sweeps in Monte Carlo simulation with the fuzzy 
X start. The dashed (dotted) line represents the 
classical value for the fuzzy x (the fuzzy S^). 



algorithm taking = Xj^'"' ' as the 
initial configuration. For the first 97,000 sweeps all the eigenvalues of Qi and Q2 are fluc- 
tuating around the classical value {Ri)^ = (^2)^ — 4.5 for the symmetric fuzzy x S^, but 
after a while the eigenvalues of Qi and Q2 arc split into the classical values {Ri)'^ — 45 and 
(i?2)^ — for the fuzzy S^, which indicates that the configuration has become the single 
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{k = 1) fuzzy S^. In this particular event we observe that the fuzzy is formed in the 
123-directions, but the formation of the fuzzy in the 456-directions should be observed 
with equal probability due to the Z2 symmetry for exchanging the two sets of directions. 

We have performed a simulation for N = 16 with the same a = 2.6 starting again from 
the fuzzy x S^, but the decay is not observed within 10,000,000 sweeps. This suggests 
that the potential barrier between the fuzzy x and the fuzzy grows with N, and 
therefore the decay probability is suppressed at large N. 
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A. Evaluation of the free energy ciround a cleissical solution 

In this section we formulate the perturbation theory and derive a formula for the free 
energy. Let us evaluate the partition function Z = J dAe~^ around the fuzzy S^ x S^ 
solution xil'^'''^'''^ at the one-loop level. The measure for the path integral is defined by 
dA = Ul=i Uali^ dAl, where A^, = Y.a=i^ with i« being the generators of SU(A^) 
normalized by tr(t"t'') = 5ah- 

We need to fix the gauge since there are flat directions corresponding to the transfor- 
mation Afj_ Afi = g A^ , where g is an element of the coset space H = U(N)/\J{k). 
In order to remove the associated zero modes, we introduce the gauge fixing term and the 
corresponding ghost term 

5g.f. = -ytr[X^,A^]2, (A.l) 
S^ost = -Ntv{[X^,c][A^,c]) , (A.2) 

where c and c are the ghost and anti-ghost fields, respectively. 

We perform the integration over A^ perturbatively by decomposing it as A^ = X/^+A^, 
where = x^^'^'^''^\ The partition function can be written as 

Z = Yo\{H)M j dAdcdce-^^"^-' , (A.3) 

where the total action -Stotai = -S" -|- Sg.f. -|- Sghost is given by 

5total = S[X] + Sun + 5int , (A.4) 

'S'kin = ^ ^ tr (^A^[Xx,[Xx,A^]]) + NtT(c[X^,[X^,c]]) , (A.5) 
+ ^ ^ a /^^p AT tr (i^i^ip) +iVtr(c[X^,[i^,c]]) . (A.6) 
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The volume of the coset space H, vol{H) = vol(U(A''))/vol(U(fc)), can be obtained by the 

^ 2 2 

formula vol(U(p)) = (p i) 1...0 1 ; a-nd the normalization factor = (27riV)~2(^ ~^ ) can be 
obtained by following the usual gauge fixing procedure as in ref. [68]. The prefactors vol(H) 
and M in (A. 3) are omitted in refs. [64,65] since they are irrelevant for the discussions in 
those papers, but we need to keep them for the comparison with the free energy in the 
Yang-Mills phase discussed in section 5. 

We calculate the free energy W = — log Z as a perturbative expansion W = Yl'j^o > 
where Wj = 0{a^^^^^^). The classical part Wq is given by eq. (2.6). Note that the kinetic 
term (A. 5) can be written as 

5kin = iVtr|^4(pA)'4 + c(pA)'c| , (A.7) 
where we have introduced the operator 

T'^mIs? [x^,M] , (A.8) 
which acts on the space of N x N traceless matrices. The one-loop term is obtained as 

Wi = 2 Tr' log {N{Vf,f} - log{vol(i/) AA} , (A.9) 

where the symbol Tr' denotes the trace in the space of TV x TV matrices omitting the zero 
modes^. 

B. Perturbative calculations for the fuzzy x 



In this section we consider the symmetric fuzzy S x S taking ni = n2 = n{= y ^) for 
the solution = xjT''''^'^^ 

B.l One- loop calculation of the free energy 

In order to solve the eigenvalue problem of the operator {V\)^, let us introduce 

= (^^.m,0y.....)®e(<^'^ (B.l) 

where Yim represents n-dimensional matrix spherical harmonics, and e^"'*) is a kx k matrix 
whose (a, b) element is 1 and all the other elements arc zero. The indices / and m represent 
the double indices (/i,^2) and (mi,m2), respectively. The matrices Y^^^''^^ form a complete 
basis of N X N matrices, and they have the properties 



''^^ ^i'm''' ^) = ^ll'^mm'^aa'Sbb' , (B.2) 



^Strictly speaking, we have to treat the zero modes for fc ^ 1 more carefully. See ref. [64] for a discussion 
on this point. This complication, however, does not affect our conclusion concerning the large- iV limit with 
fixed k. 
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Acting the operator (T-'a)^ on Y^^'^\ we find 

ir.fyt'^ = a' {hih + 1) + hih + 1)} , (B.4) 

which means that 5^^''*'* are the eigenstates of ("Pa)^- Thus the first term in eq. (A. 9) can 
be obtained as 

n-l 

2rr'log{Ar(P^)2} = 2^ J2' (2^1 + 1)(2^2 + l)log [Na^ {hih + I) + hih + 1)}] 

h,l2=0 

log-p--3 + 41ogaj , (B.5) 

where the large- hmit is taken in the second line with fixed k, and the symbol ^ ' implies 
that the zero modes = ^2 = are excluded. The second term in eq. (A. 9) can be obtained 
as 

- log I vol (If )A^} ~ (logN - ^ . (B.6) 
Adding the two terms, we obtain the one-loop contribution Wi as 

T^i = Ar2(^log^-^+41oga^ . (B.7) 
B.2 One-loop calculation of various observables 

In order to calculate various observables in perturbation theory, we need the propagators 
for An and the ghosts, which are given as 

n— 1 h I2 

;(4).,(i.k)^ = E' E E 

oh h,l2=0 mi=~li m2=—l2 

= E E' E E 

06 /i,^2=0 mi=—li m2=—l2 







(-ir.(-ir' fyt^>)Jy^„, (B.9, 



where the symbol ( • )o refers to the expectation value calculated using the kinetic term 
•Skin in (A.5) only. 

We also need to obtain the tadpole (^^), which can be expressed as 

(4) = cX^ (B.IO) 

with some coefficient c due to the SO(3)xSO(3)x Z2 symmetry. Since 

^3 3 

-tr^(X^(4)) =ctr^(X^X^) = |a2(„2_-L) , (B.ll) 

/U=l fj,= l 
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we can obtain the coefficient c by evaluating the left most term in (B.ll). At the leading 
order in we obtain 



3 6 



Itr^ (x^(4)i_ioop) = (tr^ Yl (^/^4)tr {[A,, A,][X,, A,]) \ 

/i=l \ //=1 I^,P=l / Q 

'(^""Yl Yl {Xn^tx)^^{^iaUpaA^ApA„ 

\ /i=l l/=l / Q 

in) 

Using the fact that can be written as a linear combination of ll=i,m) we can evaluate 
(B.12) as 

irS^(xiA\ \- 2 2 (2/i + l)(2/2 + l)/i(/i + l) 

(^m(^m)i-IoopJ - -^fc 2^ /l(/l + l) + /2G2 + l) ■ ^^-'^^ 

Prom (B.ll) we obtain 

^ 8fc^ (2/i + l)(2/2 + l)/i(/i + l) 

7V2a4(„2_i) ;^(Z^ + l) + l2(/2 + l) ■ ^ ■ ^ 

tl,t2=U 

Using the propagator and the tadpole obtained above, we can evaluate various observ- 
ables at the one-loop level. For instance, the "extent of space-time" is given by 



n-l 



, 6A-^(2/i + l)(2/2 + l) 

^2^ N^a''{h{h + l)+h{l2 + l)}' ^" ' 

(1,(2— U 



At large A'" with fixed a = aN i , we obtain 



The expectation value (S') can be calculated in a similar manner, but it is much easier 
to calculate it in the following way. Let us define a rescaled action 

S{X,a) = XNtv (^~[Ai„A,]'' + '^iafi,,pA^A,Ap^ (B.17) 

and the corresponding free energy 
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which is related to the original free energy W = W{l,a) through 



W{X, a) = - (AT^ - 1) log A + W{1, A3a) . (B.19) 



Then we obtain the expectation value {S) as 



{S) 1 dW{X,a) 



m dX 



A=l 



B.3 Critical point 

As we see in section 3, Monte Carlo simulations show that the fuzzy x decays into 
the pure Yang-Mills vacuum at some critical point. We can reproduce the critical point by 
perturbative calculations as follows. Let us consider the one-loop effective action for the 
rescaled fuzzy x configuration 

^ ^ 1/3(4"^ » In) for /x= 1,2,3, 

^ /3(l„®4"^)®lfc for /x = 4,5,6. 



which is given at large N as 

^r,s.xS^(/3)- Q/3^-^a/33)^ + 41og^ + logl^-^, (B.23) 

where P = I3N^ . The derivation is similar to the calculation of the free energy described 
in section B.l. (Note that one-particle reducible diagrams do not appear at the one-loop 
level.) In fact the effective action is one-loop exact in the large- A'^ limit as can be shown 
by a power counting argument [47]. The effective action has extrema at (3 satisfying 

From this we find that the effective action Fj^ 32x32 has a local minimum if and only if 
a > CKcr ^ )^ where the critical point ctcr^ ^ is given by 

-{fcS2xS2) ^4 ^^^^^1 ^ 2.481613- •• X A;3 . (B.25) 
B.4 All order results 

Since the free energy and the effective action are related to each other by the Legendre 
transformation, we can obtain the free energy by evaluating the effective action at its 
extremum. Then exploiting the fact that the effective action is one-loop exact, we can 
evaluate the free energy to all orders in perturbation theory. This method is proposed in 
ref. [34], and it is applied to other models successfully [65,67]. 

Above the critical point a > Ocr ^ \ the value of P that gives the local minimum of 
the effective action can be obtained by solving eq. (B.24) as 



^ = /(«) = f (1 + V1TI+ j2-<5+^===j , (B.26) 
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where 



At large a, the solution (B.26) is expanded as 

- , / 4/c 48A;2 960k^ \ , 
/3 = / a = a 1 - ^ - ^3 . B.28 

Plugging this solution into (B.23), we obtain the free energy to all orders as 
1 rx. a"^ _ , 16A^^ 15 8k 224^ 3520^^ 

— -:r Wh^2^o2 ~ + 41oga + log —-Z J s r;^ • • • . (B.29) 

Using (B.20), we obtain the all order result for the expectation value {S) as 
1 5 8A: 448^^ 3520^^ 

We can also calculate other observables to all orders as in ref. [65]. For instance, let 
us calculate (;^tr(yl^)^). Since the one- loop contribution comes totally from the tadpole 
diagrams, which are one-particle reducible (The second term in eq. (B.15) is subleading 
in N.), we may obtain the all order result by simply replacing a by the quantum solution 
/(a) in the classical result as 

1 / 1 , , .nX 1 ,2 4 40k 768A:2 18304/c3 



^\N / ~ 2k-'^ ' 2k a2 ^6 c^io 

C. Perturbative calculations for the fuzzy 

In this section we briefly describe the perturbative calculations for the fuzzy taking 
ni = n(= ^) and n2 = 1 for the solution = x^"^'"^''^^ 

C.l One- loop calculations 

The one-loop free energy is given at large AT by 

-^W,s^ = -^(n'-l) + 2k'J2(^l + l)log[NaH{l + l)]-\og{vol{H)M] 

1=1 

where a = aN^ . The tadpole is given by 



(A). , - i for//-l,2,3, 
Wi-ioop-<^Q for /x = 4, 5, 6. ^^-^^ 
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The observables can be evaluated as 



n\n 



1— loop 



trQi^ ^ = ^{tr(X^X^) + 2tr (x^(4)i-ioop) + (tr(i^)')o} 

2ij_. 2 ^^ 4fc^(n'-l) 6 "^ 21 + 1 
" 1 4A^^" ' N^a^ ^ Nv?a^ ^1(1 + 1 



1 ,^2 -iN 4A: ^(n^ - 1) 6 2^ + 1 

a2 4 

^ ^ ^ trQ2^ ^ , (C.4) 

^,4 



^(5)i-ioop = -^(n^ - 1) + ^ + 4^{6(-^' - 1) + 



24fc2 2 



C.2 Critical point 



The critical point can be obtained by considering the effective action for the rescaled fuzzy 
configuration 

^ ^ f /3(4"^ Ifc) for// = 1,2, 3, .^g. 
^ \q for /X = 4, 5, 6 . ^ ■ ^ 

The effective action is obtained in the large- A'" limit as 



where (3 = (3\ N . The critical point is obtained as 

4 



^ 243^^ = 2.9511518- •• x Vk . (C.8) 



C.3 All order results 

Above the critical point a > ctcv , the effective action has a local minimum at 



^ = 5(5) = 1(1 + ^1+1+ W2-£ + -^j, (C.9) 



£ = a-3(256A;2)3 < 



1 , 1 

3 



The all order result for the free energy can be obtained by substituting /? = g{d) in (C.7). 

1 ~ , 11 16A;2 896/0^^ 28160^^ , 
]^^.s^--^+41oga + log^-^-^-^-^^X^ -^^-^l) 
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The expectation values of observables can be calculated at large N as 



1 a"^ 5 16/c2 1792/c^ 28160/c^ 

n\ 4 80fc^ 3072fc^ 146432fc6 

^(^trg,)^0. (C.14) 

D. Free energy in the Yang-Mills phase 

In ref. [71] the free energy for the pure Yang-Mills model (a = 0) is calculated as 

^i^-i^|/z. + ^ln(^/fiv)-ilnvr| , (D.l) 

where D is the number of bosonic matrices, which is = 6 in the present case, and the 
third term is added to adjust the normalization of the measure to the one used in the 
present paper. The "free energy density" fjj can be calculated analytically by the 1/D 
expansion [49], and one obtains /oo = — l/4atD = oo. For D = 6 the gaussian expansion 
method [71] gives /g = —0.33 (the convergence is clear up to order 7). Therefore the free 
energy in the Yang-Mills phase of our model is given at large N as 

-^WVM^31ogA^ + 6 {-0.33 + log(3i7r-^)} , (D.2) 

if a is sent to zero in the large-iV limit. Finite a corrections should be 0(0;^) due to the 
parity symmetry of the pure Yang-Mills model. 
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